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Abstract 

This paper is concerned with the discrete spectra of Schrodinger operators 
H = —A + V, where V(r) is an attractive potential in N spatial dimensions. 
Two principal results are reported for the bottom of the spectrum of H in each 
angular-momentum subspace Tig : (i) an optimized lower bound when the poten- 
tial is a sum of terms V(r) = V^\r) + V^ 2 \r) , and the bottoms of the spectra of 
_A+l/( 1 )(r) and -A+V (2 \r) in Til are known, and (ii) a generalized comparison 
theorem which predicts spectral ordering when the graphs of the comparison poten- 
tials V (1) (r) and V {2) (r) intersect in a controlled way. Pure power-law potentials 
are studied and an application of the results to the Coulomb-plus-linear potential 
V(r) = —a/r + br is presented in detail: for this problem an earlier formula for 
energy bounds is sharpened and generalized to iV dimensions. 
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1. Introduction 

This paper has two principal aspects: the potential-sum approximation, and 
the generalization of the comparison theorem of quantum mechanics to cases where 
the comparison potentials intersect. We study spherically-symmetric problems in N 
spatial dimensions. There is much interest in problems posed in arbitrary dimension 
N [1-9], rather than specifically, say, for N = 1, or N = 3. References [1] and [4] 
are useful for technical results such as the form of the Laplacian in iV-dimensional 
spherical coordinates; the other papers are concerned with solving problems such as 
the hydrogen atom [3, 7] and the linear, harmonic-oscillator, hydrogen atom, and 
Morse potentials [8] in higher dimensions than N = 3. The geometrical methods 
we use in this paper are independent of dimension, which can usually be carried 
as a free parameter N. We consider examples with Hamiltonians of the form, 
H = —A + v sgn(q)r q , or with sums of such potential terms. We suppose that the 
Hamiltonian operators H = — A + V(r), r = ||r||, have domains V(H) C L 2 (R N ), 
they are bounded below, essentially self adjoint, and have at least one discrete 
eigenvalue at the bottom of the spectrum. Because the potentials are spherically 
symmetric, the discrete eigenvalues E n £ can be labelled by two quantum numbers, 
the total angular momentum t = 0, 1, 2, ... , and a 'radial' quantum number, n = 
1, 2, 3, . . . , which counts the eigenvalues in each angular-momentum subspace. Since 
the discrete spectrum may be characterized variationally [10,12], the elementary 
comparison theorem < E^} < E$ immediately follows. The 

generalization we shall study (in Section 3) involves comparison potentials whose 
graphs 'cross over' in such a way that spectral ordering is still guaranteed. 

Before we study the generalized comparison theorem, we shall need some es- 
tablished results concerning 'kinetic potentials' [13] and 'envelope theory' [14, 15]. 
In order to fix ideas and simplify the presentation, let us suppose that E is a dis- 
crete eigenvalue at the bottom of the spectrum of H = —A + V in N dimensions. 
It follows that E = M(ip,Hip) where ip G V(H), and ||^|| = 1. We perform 
the total minimization in two stages: first we constrain the process by fixing the 
mean kinetic energy (ip, —Aip) = s, and then we minimize over s > 0. The mean 
potential-energy function under the constraint is called the 'kinetic potential' V(s) 
associated with the potential V(r). Thus we define 



V(s) 




{s + V(s)}. 



(1.1) 
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The variational definition of the kinetic potentials implies that (i) cV(s) = cV(s), 
and (ii) < V {2) (s) E^ < E^ . Kinetic potentials can be defined [13] for 

higher eigenvalues and they can then be reconstructed from 'energy trajectories', 
the functions which describe how the eigenvalues vary with the coupling parameter 
v > 0. We have in general for coupling 

H = -A + vf(r) -> Eru = F n£ (v) (1.2a) 

and 

a = F n£ (v) - vF^(v), J n£ (s) = FUv), (1-26) 

where E n£ is the n th eigenvalue in the angular-momentum space labelled by £, 
and F n£ (v) describes how this eigenvalue depends on the coupling v > 0; the 
corresponding kinetic potentials may then be defined by (1.2b). The relationship 
F(v) <-> f(s) is essentially a Legendre transformation [16]: for the ground state (or 
the bottom En of each angular-momentum subspace) F(y) is concave [11, 12, 
13] and consequently f(s) is convex: it follows [17] immediately from (1.2b) that 
F"(y)f"(s) = —v~ 3 < whenever F"(v) ^ 0; thus in general F(y) and f(s) 
have opposite convexities almost everywhere. For the important class of examples 
H = —A + v sgn(g)r g , corresponding to pure powers q > —2, we know that 
F(v) is concave for every (discrete) eigenvalue since, by scaling arguments, we have 
F n£ (v) = F n£ (l)v^ , and sgn(F n ^(l)) = sgn(g). 

The main purpose for this two-step reformulation of 'min-max' is that certain 
spectral approximations are very effectively developed in terms of kinetic potentials. 
We shall consider first the 'envelope approximation', which in its most succinct form 
can be summarized as follows 

f(r)=g(h(r)) U{s) » g(K~ £ (s)). (1.3) 

Here f(r) is a smooth transformation of a 'base' potential h(r). We suppose that 
the transformation g is monotone increasing and, if it also has definite convexity, 
the following important conclusions may be drawn: if g is concave, we get an 
upper bound ~=<; and, if g is convex, we obtain a lower bound ^=> . These 
results may also be derived by the use families of upper and lower 'tangential' 
potentials [18]. In Section 2 we shall apply this result to study the Coulomb-plus- 
linear potential V(r) = —1/r + Xr which is clearly at once a convex transformation 
of the Hydrogenic potential h(r) = — 1/r and a concave transformation of the linear 
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potential h(r) = r. We shall show that we are also able to express both the upper 
and lower bounds for the entire discrete spectrum in the form of explicit rational 
functions A = X(E n i). 

The base potentials used for the Coulomb-plus-linear potential are both pure 
powers. Thus we shall need to use the corresponding base kinetic potentials. In fact 
we have shown in general [14] that 

-A + sgn(g)r« =► P^ = nunj(^M) +sgn(^|, (1.4) 

where, for example, P^(-l) = (n+i+N/2-3/2), and P$(2) = (2n+£+N/2-2). 
These P-numbers and the underlying eigenvalues E^ e satisfy the relation E^ e = 
En0~ 2£ '■ this is generally true for central potentials and is the content of Theorem 2, 
which we prove in Section 4. Numerical values for P^(l) are given in Table 1 for 
N = 2, 12 . It is interesting that the case q = corresponds exactly to the 
ln(r) potential [15]. The expression in (1.4) is derived by a change of variable 
s — > (P N /r) 2 in the kinetic-potential formalism. The application to the Coulomb- 
plus-linear potential is not our only interest in these P^-numbers. They provide 
through (1.4) a nice representation for the pure-power eigenvalues since the P N - 
numbers vary smoothly with q through q = whereas the eigenvalues themselves 
do not [14]. We have proved [14] that P^{q) are monotone increasing in q. This 
result was obtained by using envelope theory: we considered one power q as a 
smooth transformation of another p , and then took the limit p — > q. 

In Section 3 we prove Theorem 1 which provides a lower bound for the bottom 
of the spectrum in each angular momentum subspace using the sum approximation. 
In Section 4 we prove Theorem 2, which establishes the invariance of the eigenvalues 
with respect to changes in £ and N that leave the sum N+2£ invariant. This allows 
us to restrict our considerations to the ground state in sufficiently high dimension 
N. We reformulate the refined comparison theorem (Theorem 3 of Ref. [19]) 
which becomes Theorem 3 here. We first prove the monotonicity of the ground- 
state wave function in N dimensions; then we prove Theorem 4, which extends 
Theorem 3 to N > 2 dimensions. Finally we prove Theorems 5, 6, and 7 which 
provide simple explicit sufficient conditions for the application of Theorem 4 under 
a variety of crossing schemes. In Section 4 we apply Theorem 5 to sharpen the 
envelope bounds already found in Section 2 for the bottom of the spectrum E of 
H when V is the Coulomb-plus-linear potential V(r) = —a/r + br. 
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2. Coulomb-plus-linear potential: an eigenvalue formula 

The Coulomb-plus-linear potential V(r) = —a/r + br is of interest in physics 
because it serves as a nonrelativistic model for the principal part of the quark-quark 
interaction. First, we will use the envelope method to derive a simple formula for 
upper and lower bounds for all the eigenvalues E n g 7 n=l,2,3,..., £ = 0,1, 2,... 
Because the linear potential, rather than the harmonic oscillator, is used as a basis 
for the upper bound, the new bounds are sharper than those of the earlier paper [20] . 

If we denote the eigenvalues of H = —uiA — a/r + (3r by E{uj,a,f3) and 
consider a scale of change of the form r' = r/a , and if we further choose a = a/ui , 
then it is easy to show that 

Buj 2 

E(u, a, (3) = a^-^fl, 1, A), A = ——. (2.1) 

a 6 

Thus it is sufficient to study the special case H = — A — 1/r + Ar. 

We need a solvable model which we can use as an envelope basis. The natural 
bases to use in the present context are the hydrogenic and linear potentials 

h(r) = sgn((/)r 9 , where q = —1, 1. (2.2) 

The spectrum generated by the potential h(r) is represented precisely by means of 
the semi-classical expression (1.1) as follows: 

£„t{v) = min{s + vh n£ (s)}, (2.3) 

s>0 

where the 'kinetic potentials' h n i(s) associated with the power-law potentials (1.1) 
are given [15] by 

g + 2 
2 

s~ q / 2 , (2.4) 

and S^g is the eigenvalue of —A + sgn(q)r g in iV dimensions, that is to say, 
corresponding to the pure-power potential with coupling 1. If we use the potential 
h(r) = — i as an envelope basis, then V(r) = + Xr = g{~) implies g is 
convex. And if we use the linear potential h(r) = r as an envelope basis, then g is 
concave. A weaker upper bound is provided by the harmonic oscillator h{r) = r 2 , 
for which, again g(h) is convex. 

For the power-law potentials h(r) = sgn(q)r q we can simplify (2.3) by changing 
the minimization variable s to r defined in each case by the equation h n i(s) = h(r) 



hni(s) 



2 + q 
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so that g(h(r)) = f(r) = — - + Xr . The minimization on the other hand, which 
yields eigenvalue approximations for the Hamiltonian H = —uA + f(r) (u> > 0), 
can be expressed in the form 



E. 



N 



mm 

r>0 



h Xr 

r 



(2.5) 



where 



pN 



E. 



(<?) 

nl 



2 + q 
2q 



2 + q 



2 + q 



(2.6) 



We obtain a lower bound with P^(— 1) = (n + £ + N/2 — 3/2) and the harmonic- 
oscillator upper bound (of Ref.[20]) with P$(2) = 2n + £ + N/2 - 2, and a sharper 
upper bound with P^(l); the P,fo(l)-numbers are provided in Table 1 for N = 



2, 12. This table allows £ > since P* e = P 



>N+2e 

nO 



it is clear that E^ £ (—l) and 
E^ f (2), and the corresponding P-numbers, are invariant with respect to changes in 
£ and N which preserve the sum 2£ + N; this symmetry is also true for £^(1), 
indeed for all eigenvalues generated by a central potential. This property is the 
content of Theorem 2, which we state and prove is Section 4. We thus obtain the 
following energy bounds 



mm 

r>0 



- + Xr) > < E n g < min 

r I r>0 



+ Xr 



(2.7) 



for n = 1, 2, 3, £ = 0, 1, 2, .... Consequently, the energy bounds are given by the 
parametric equations 

E„.f = - — + ^ (2.8a) 



2vt 2 

1 = ^ + ^, t = rP»{ q ) , = -1,1, 



(2.8b) 



wherein the lower and upper bounds take the values v = fx = P^(— 1) and v = 
[i = P^(l) respectively. It is interesting that we can actually solve Eqs. (2.8a) and 
(2.8b) to obtain A as an explicit function of E = E^ £ ; the result namely is 



A = 



2{vEy 



vE< 



(l + 3v 2 E)? -l]} 



(1 + 3z/ 2 £)5 - 1 



(2.9) 
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with E > —4^2 (corresponding to A = for the pure hydrogenic spectrum). We 
emphasis that these bounds are valid for all the discrete eigenvalues in arbitrary 
dimension N > 2. The bounds are weak for n > 1 , but at the bottom of each 
angular momentum subspace n = 1 they are sharp and improve with increasing 
£, N, and A. The lower bound for the bottom of each angular-momentum subspace 
(n = 1) can be improved by use of the 'sum approximation' ( [21] and Section 3 
below) in which v = P%(1) (Table 1) and fj, = P%(-1) = (n + £ + N/2 - 3/2) . In 
Figure 1 we exhibit these bounds for n = 1 , N = 3, and £ = 0, 1, 2, 3. 

3. The sum approximation: lower bounds 

We now consider potential which are sums of terms. Since further generaliza- 
tions easily follow, we first look at the problem of the sum of only two potential 
terms. We assume that each potential vhS l \r) alone, when added to the kinetic- 
energy operator —A, has a discrete eigenvalue E at the bottom of the spectrum 
for sufficiently large coupling v. We note that the proof is unchanged if we restrict 
the problem to a given angular-momentum subspace labelled by £; our claim then 
concerns the bottom of the spectrum of H in such a subspace; in the more general 
case, all the kinetic potentials would be labelled by £. We express our result in 
terms of kinetic potentials and prove (for the case £ = ) the following: 



Theorem 1 If E is the bottom of the spectrum of the Hamiltonian H = —A + V, 
and the potential V is the sum V(r) = hM(r) + /i (2) (r), then it follows that the 
sum of the component kinetic potentials yields a lower bound to V, that is to say 

V(s)>h (1 \s) + h (2 \s). (3.1) 

We shall now prove this theorem, which is in effect an optimized Weyl lower 
bound [22-24]. From the definition (1.1) of kinetic potentials we have 

V(s) = inf (ib, V%b) = inf (V, (h {1) + h i2) ) ip] . 

(V,V)=i W>,V0=i 

(tP,K4>)=s {il>,Kip)=s 

But the latter minimum mean-value is clearly bounded below by the sum of the 
separate minima. Thus we have 

V(s)> inf (ip,h (1) ^)+ inf (V,>i (2 V) =h {1 \s) + h (2 \s), 
(V,V)=i W>,V0=i 

(ip,-Aip)=s (ip,-Aip)=s 
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which inequality establishes the theorem. □ 

Another approach, which would eventually yield an alternative proof of the 
theorem, exhibits the relationship between Theorem 1 and the classical Weyl lower 
bound [24-26] for the eigenvalues of the sum of two operators. Let us suppose that 
\1/ is the exact normalized lowest eigenfunction of H = — A + V, so that H^> = E^>. 
If the positive real parameter w satisfies < w < 1, then E = (—A + V)^) 
may be written as follows: 

E = w(v, (-A + 1 h {1) (r)j ^ + (1 - w) (v, (-A + -^-^ h i2) (r)^j ^ 

> w inf ( if), ( -A + - hW(r) J V J 

ipET>(H) V. V w J J 

(V,^)=i 

+ (1-w) inf (^(-A+—}—hW(r))ifl). 

That is to say, in terms of component kinetic potentials, we arrive at Weyl's in- 
equality for the lowest eigenvalue E of the operator sum H, where 

H = -wA + h {1) + -(l-w)A + h (2 \ 

and we conclude 

E > w min { s H }i \s) \ + (1 — w) min { s H }i Vs) 1 . 

~~ s>0 { W J s>0 [ 1-W J 

Since w is an essentially free parameter in the last expression, we may optimize 
the Weyl lower bound with respect to the choice of w : this forces the individual 
values of s at the minima, {si(w), S2(w)}, to be related. More specifically we find 
from the individual minimizations over s, 

E > £(w) = w Sl {w) + (1 -w)s 2 (w) + h (1 \s 1 (w)) +h (2 \s 2 (w)), 



where 



dh {1) dh {2) 

w = --fe M w ))> and l ~ w = ~~d^ ( s z( w ))- 



The critical condition S'(w) = for the subsequent maximization over w then 
yields si(w) = S2(w). Thus the best lower energy bound is given by 



E > min|s + /i (1) (s) + /i (2) (s)}. 
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The kinetic-potential inequality of Theorem 1 leads, of course, to the same energy 
lower bound: the optimization just performed above is therefore seen to be auto- 
matically 'built in' by the kinetic-potential formalism. 

It follows immediately from the above kinetic-potential comparison theorem 
and coupling-parameter absorption that a lower bound to the lowest energy E of 
the Hamiltonian H = —A + ^\ Cih^\r), {ci > 0}, is provided by the formula 

E > mm js + ^c^ ( %)j . (3.2) 

Similarly we can extend this result to 'continuous sums' such as V(r) = 
£c(t)hU(r)dt. 

Meanwhile, since the proof is identical, the bound is valid for the bottom of 
each angular-momentum subspace. Thus, more generally, the fundamental inequal- 
ity becomes 

V ie (s)>h^(s) + hf £ \s), £ = 0,1,2,.... (3.3) 



4. Generalized Comparison Theorems 

The proof of our generalized comparison theorem (Theorem 4) depends on 
monotone behaviour of the wave function induced by the assumed monotonicity of 
the potential. We are able to establish this monotonicity for the lowest eigenfunction 
in arbitrary many spatial dimensions N > 1. We shall then be able to apply our 
eigenvlaue results to the case £ > and n = 1 because of Theorem 2 which claims 
that E^ £ = E^Q~ 2t ; this general result is then employed in the special case n = 1. 
Theorem 2 

Suppose that H = —A + V(r), where V(r) is a central potential in N > 2 dimen- 
sions, has a discrete eigenvlaue E^ £ with n radial nodes in the angular-momentum 
subspace labelled by £, then E^ £ = E^ 2i . 

Proof: We suppose that ip is the eigenfunction corresponding to E^ e . We express 
—A in spherical coordinates [1-9] and write the radial eigenequation explicitly as 

-V» - ^^W) + * ( * + r T~% (r) + V{r)^{r) = E^(r). 

If we now define the reduced radial function u(r) e L 2 (R + ) by ip(r) = 
u(r)r- (Ar-1)/2 , r > 0, and u(0) = 0, we obtain 
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-u"(r) + 



(JV - 1} (N 2 3) + *V + N - 2) 



+ V(r) 



u(r) = £>(r) 



(4.1) 



If we consider the spherically-symmetric potential V(r) in M dimensions such that 
(M - 1)(M - 3)/4 = + JV - 2) + (N - 1)(N - 3)/4 , we find that M = 2£ + N. 
The eigenequation (4.1) then may be written equivalently 



-u"(r) + 



(M-l) (M-3) 
2 2 



+ V(r) 



(4.2) 



It therefore follows immediately that = E™ = E% +N . □ 
For the purpose of our comparison theory we may now consider the special 
case n = 1, £ = in arbitrary N > 1 spatial dimensions: the energy results 
which we derive will then be applicable to the family of equivalent problems in N' 
spatial dimension with n = 1, £ > 0, and N = N' + 2£. In order to prove an 
appropriate extension of the comparison theorem in iV dimensions, we shall first 
need to establish an elementary monotonicity property for the ground-state ip. We 
prove the following: 

lemma 

Suppose ip = ip(r), r = \\r\\, r E R N , satisfies Schro dinger's equation: 



Hi(;(r) = (-A + V(r))^(r) = Eip(r), 



(4.3) 



where V(r) is a central potential which is monotone increasing for r > 0. Suppose 
that E is a discrete eigenvalue at the bottom of the spectrum of the operator H = 
—A + V, defined on some suitable domain T>(H) in L 2 (R N ). Suppose that ip(r) 
has no nodes, so that, without loss of generality, we can assume ip{r) > 0, r > 0. 
Then ^'(r) < 0, r > 0. 

Proof: The proof for the case iV = 1 is given in Ref. [19], Eq.(2.2). Henceforth we 
shall now assume N >2. If we express —A in spherical coordinates in iV spatial 
dimensions, then we have 

-Aip + Vip = Eip 
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We now multiply by t both sides and integrate with respect to t , to obtain 

<//( r ) = (l/r^" 1 ) /V(t) - E\il)(t)t N - 1 dt. 
Jo 

Since y is monotone increasing, it follows that there is one point f > satisfying 
V(r) = E. First, we prove that s(r) = /- [V(t) — E]ip(t)t N ~ 1 dt is monotone in- 
creasing and bounded. For t > f, [V(t) - E^it)^' 1 > 0, because V(t) > E 
and hence s(r) is increasing as r — > oo . If there exists ri < oo such that 
s( ri ) = - | Q f [V(t) - E\il)(t)t N - 1 dt , then J^[V(t) - E)^{t)t N - x dt > 0, r > n , and 
tfj'(r) > 0, r > T\ ; this contradicts the fact that the wave function ip(r) is positive 
and belongs to L 2 (R N ) . This means that ^[V(t) - E^^t^dt < 0, Vr > . 
Consequently ifj'(r) < 0, Vr > 0. □ 

We now consider two potentials V\{r) and V 2 (r) both of the type described 
above. We have two Schodinger equations for the respective ground-states and 
ip2 and the corresponding discrete eigenvalues E\ and E<i at the bottoms of the 
spectra. Thus we have the following pair of eigenequations 

(-A + Vi(r)tyi(r) = EtMr) (4-4) 

(-A + y 2 (r))^ 2 (r) = E 1 Mr) (4.5) 

The radial wave functions in the present paper satisfy the normalization condition 
J oo ^|(r)r Ar - 1 rfr < oo, i = 1,2. With this notation, and iV = 3, Theorem 3 of 
Ref. [19] becomes 
Theorem 3 

k(r) = [ (Vi(£) - V 2 {t))^i{t)t 2 dt < 0, Vr > 0, i = 1 or 2 =>- £ x < £ 2 - (4.6) 
Jo 

We shall now generalize this theorem to general dimension N > 1 . We first establish 

a fundamental comparison formula (Eq.(4.7)) below. 

By multiplying (4.4) by ip 2 and (4.5) by ifti, and subtracting, we find 

- V^A^i + [Vi - V 2 ]ipiip 2 = [Ex - E 2 }ipiip2 
Integrating over R N and using the following identity, 
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v.(v>iVV> 2 ) = W>i.v-i/>2 + ^iV 2 -02 

we find that, 



/ V#iVV> 2 - ih^i]d N r + [ [Vi - V 2 ]^ 2 d N r = [E x - E 2 ] [ ipnp 2 d N r 
Jr n Jr n Jr n 

Now by Gauss's theorem [27] we find that the first term becomes a surface integral 
which vanishes because ipi G L 2 (R N ). In the remaining integrals the angular factors 
yield 2n N / 2 /T(N/2). Hence we find 



2vV 



[Vi(r)-V 2 (r)]Mr)Mr)r N - 1 dr=^^[E 1 -E 2 ] J ^(r)^^)^" 1 *, 



r(iv/2) J 

which implies, 



f oo 

S 



/ [Vi - V 2 ]ipiip2r N - 1 dr = [E x - E 2 ] / ^ 2 r N ~ l dr (4.7) 
Jo Jo 

Now we may state our generalization of Theorem 3 to N dimensions: 

Theorem 4 

k(r) = [ (VI (t) - V 2 (t))i; i (t)t N - 1 dt < 0, Vr > 0, i = 1 or 2 ^ £?j < £? 2 . (4.8) 
Jo 

Proof: For definiteness we assume that i = 1; the proof is just the same with the 
other choice. We study the integral s on the left side of (4.7) . Integrating by parts 
we find that 

/•OO 

s = [fc(r)V2(r)]g° - / k(r)iP' 2 (r)r N - 1 dr (4.9) 

Since fc(0) = ip 2 (oo) = 0, the first term vanishes, and s is therefore equal to the 
negative of the integral of the right side of (4.9) . But the integrand of this integral 
is positive because k(r) < 0, by hypothesis, and we know that ip' 2 (r) < by the 
above lemma. This proves that s < 0. Consequently, by (4.7), we obtain E\<E 2 . 
□ 

It may be difficult to apply Theorem 4 in practice. Thus it would be helpful to 
establish some simpler sufficient conditions, depending on the number and nature 
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of the crossings over of the two comparison potentials. We treat three useful cases: 
Theorem 5, one potential crossing, with use of the wave function; Theorem 6, two 
crossings and the use of the wave function; Theorem 7, two crossings and no wave 
function used. In these Theorems we shall assume that the integrals f Q (Vi(r) — 
V 2 (r))ifji(r)r N ~ 1 dr, % = 1,2, exists for the given problem, even though we use at 
most one wave function factor. 

Theorem 5. If the potentials V\{r) and V 2 (r) cross exactly once for r > at 
r = r\ , with, 

(i) Vi(r) < V 2 (r) (0 < r < n) and 

(™) Jo°°[W) - V 2 (t)]Mt)t N - 1 dt < 0, i = 1 or 2, 

then 

k(r) = [ [Vi(£) - V 2 (t)]t N - 1 ^ i (r)dt < 0, V r > 0, % = 1 or 2, (4.10) 
Jo 

from which E\ < E 2 follows, by Theorem 4- 

Remark: The best bound is obtained with the equality in hypothesis (ii). 

Proof of Theorem 5: We choose i = 1 : the proof is identical for % = 1 or 2. 
First, we show that s(r) = [Vi(t) — V 2 (t)]ipi(t)t N ~ 1 dt is monotone increasing. 
For t > ri, s'(r) = [Vi(r) - V 2 {r)]tl) 1 (r)r N ~ 1 > 0, because Vi(r) > V 2 (r); hence 
s(r) is increasing on (ri,oo) . Moreover, (ii) implies that the maximum value of 
s(r) is reached at r = oo; i.e s(r) < s(oo) we have therefore 

/>oo 

/ [v 1 (t)-v 2 (t)]Mt)t N - 1 dt = 

Jo 

/ [v 1 (t)-v 2 (t)]Mt)t N - 1 dt+ m^-v^Mt^-'dtKo 

Jo Jr x 

and therefore 

lim s{r) < - f \vtit) - V 2 (t)]1> 1 (t)t N ~ 1 dt- 
Now, we have the following two cases to consider 

Case 1: for r < n , k(r) = ^[Vxit) - V^ip^t^^dt < 0, since V t (t) < V 2 {t) 
for < t < r. 
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Case 2: if r > r\ , then 

k(r) = [ 1 [V 1 (t)-V 2 (t)]if; 1 (t)t N - 1 dt+ [ [V^tyViWMt)**' 1 ^ = s(r)-s(oo) < 

JO Jr x 

Therefore, k(r) < 0,Vr > 0. □ 

Theorem 6. If the potentials V\(r) and V 2 (r) cross twice for r > at r = 

ri , r = r 2 (r x < r 2 ) witt, 

(%) Vi(r) < V 2 (r) /or < r < r x and 

(it) fl\V x {t) - V 2 {t))^ l (t)t N - 1 dt < 0, z = 1 or 2, 

t/ien, 

fc( r ) = / (^(t) - V 2 (t))i; i (t)t N - 1 dt < 0, Vr > 0, i = 1 or 2, (4.11) 
Jo 

from which E\ < E 2 follows, by Theorem 4- 

Remark: The best bound is obtained with the equality in hypothesis (ii). 

Proof of Theorem 6: k'(r) = iV\{r) — V 2 (r))r N ~ 1 ipi(r) , now k(0) = 0, k'(r) < 
0, < r < ri, implies k(r) < 0, < r < r\. Next, k{r 2 ) = 0, k'(r) > 0, r\ < 
r < r 2 , implies k(r) < 0, r\ < r < r 2 . Lastly, k(r 2 ) = 0, k'(r) < 0, r > r 2 , 
implies /c(r) < 0, r > r 2 . □ 

Theorem 7. If the potentials V\{r) and V 2 {r) cross twice for r > at r = r\, r 2 

(ri < r 2 ) 

(%) Vi(r) < V^r) /or < r < n and 

r«; ^(v^-v.m^dtKo 

then, 

k(r) = [ [Vi(£) - V^)]^)** -1 *** < 0, Vr > , i = 1 or 2, (4.12) 
/rom which E\ < E 2 follows, by Theorem 4- 

Remark: The best bound is obtained with the equality in hypothesis (ii). 

Proof of Theorem 7: We choose i = 1 : the proof is identical if i = 2. Define 
h(r) = J r (Vi(t) - V 2 (t))t N - 1 dt the proof of Theorem 3 shows that h(r) < 0, < 
r < r 2 . But 

Jo 
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= [h(t)Mt)] r o - / Ht^[(t)dt 
Jo 

= h(r)ip 1 (r)- f h(t)ip[(t)dt < 0, 
Jo 

meanwhile, k'(r) < 0, r > r 2 . Therefore, k(r) < 0, Vr > 0. □ 

5. Application to the Coulomb-plus-linear potential 

As an example, we employ the comparison theorems to improve the bounds 
obtained in Section 2 for the eigenvalues corresponding to the Coulomb-plus-linear 
potential V(r) = —a/r + br, where a and b are positive coupling parameters. 
For the upper bound we use as a comparison potential the shifted linear potential 
h(r) = —a + fir, where a and (3 > . We allow the potentials V(r) and h(r) 
to cross over exactly twice, as illustrated in Figure 2. Let A and B represent the 
absolute values of the (or of the -weighted ) between the potentials. We 
vary a and (3 so that A = B , and thereafter Theorems 5 and 6 imply E v < E h . 
For simplicity of derivation of the upper-bound formula, we will use Theorem 7 
(with no use of the wave function ifj ) . Thus we have two equations to solve in this 
case, 

a 7 n 

h br = —a + pr, 

r 

/ [ — + bt + a- f3t]t N ~ 1 dt = 0, 
Jo t 

where r = r 2 is the second crossing point. These reduce to the quadratic equations 

(J3 - b)r 2 - ar + a = 0, 

iV(iV - 1)(6 - /3)r 2 + a(N - 1)(N + l)r - aN(N + 1) = 0, 

with simultaneous solution r = • Now the best upper bound is obtained 

after minimizing with respect to r, giving 

^{-G^M^O'H (5 - 1) 

At the expense of further complication, the use of ipi (r) (the Airy function) would 
lower this upper bound. 

Similarly, to improve our lower bound, we allow the Coulomb-plus-linear poten- 
tial to intersect twice with the Hydrogenic potential h(r) = —f+/3, with the exact 
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wave function solution ip = e - ar /( Ar - 1 ) anc [ the exact energy E h = f3—a 2 / (N—l) 2 , 
where a and (3 are positive parameters. Again, let A and B represent the ab- 
solute values of the areas (or of the ift -weighted areas) between the potentials. We 
vary a and (3 so that A = B , as illustrated in Figure 3, and thereafter Theorems 5 
and 6 imply E h < E v . Subsequently, we obtain the lower bound (without ifj ) by 
solving the following three equations: 



Jo 



—a , a 

-T + bt = -J + P 



E L = mm{P-(a/(N-l)) 2 }. (5.2) 

For the case a = 1 and b = 1 , we compare in Figure 4 the upper and lower 
bounds obtained by means of the classical envelope method and by the comparison 
theorems introduced in Section 3. Generalizations to cases where there are a large 
number of potential crossings are discussed in Ref. [37]. 



6. Conclusion 

Our proof of the lower-bound for the bottom of the spectrum of the operator 
H = —A+V^(r)+V^(r), based on kinetic potentials, is more compact and direct 
than the original proof, and is valid in N dimensions; the principal steps of the 
earlier proof are repeated because they show that the final result is equivalent to an 
optimization of the classical theorem of Weyl. The generalized comparison theorem 
is proved in the present paper for all dimensions N, whereas, in its original form, 
it required two distinct theorems, for N = 1, and N = 3. Moreover, we are 
now able to apply the results to the bottom of each angular-momentum subspace 
since we have proved that this energy is identical to the lowest eigenvalue of a 
higher-dimensional problem, in R N+2£ . Meanwhile, in order to be practical, weaker 
sufficient conditions were sought which would guarantee in a simple way that the 
comparison potentials cross over so as to imply definite spectral ordering. These 
results greatly clarify the application of the generalized comparison theorem to 
specific problems. 

The Coulomb-plus-linear problem provides a convenient example on which to 
test the effectiveness of the energy bounds. At the same time it offers an opportunity 
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to sharpen an earlier energy-bound formula for this problem, and to extend its va- 
lidity to all N > 2 dimensions. For most of the parameter space of the problem, the 
energy bounds provided by this formula for the bottom of each angular- momentum 
subspace (n = 1) are accurate to a few percent and, as we have shown, they become 
sharper with increasing N or t. If the sum approximation is capriciously applied 
also to the higher discrete eigenvalues n > 1 , the resulting ad hoc approximation 
formula continues to give very accurate estimates, which, however, are no longer 
bounds. What additional conditions might guarantee bounds from such a formula 
is an interesting open question. 
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Table 1 The 'input' P-values -P^o(l) used in the general formula (1.4), for N = 
2, 3, . . . , 12 . The same data applies to £ > since, by Theorem 2, we have P^ e — 



pN+2£ 
r n0 



N 


n = 1 


n = 2 


n = 3 


n = 4 


2 


0.9348 


2.8063 


4.6249 


6.4416 


3 


1.3761 


3.1813 


4.9926 


6.8051 


4 


1.8735 


3.6657 


5.4700 


7.2783 


5 


2.3719 


4.1550 


5.9530 


7.7570 


6 


2.8709 


4.6472 


6.4398 


8.2396 


7 


3.3702 


5.1413 


6.9291 


8.7251 


8 


3.8696 


5.6367 


7.4204 


9.2129 


9 


4.3692 


6.1330 


7.9130 


9.7024 


10 


4.8689 


6.6299 


8.4068 


10.1932 


11 


5.3686 


7.1274 


8.9053 


10.7453 


12 


5.8684 


7.6253 


9.4045 


11.2744 
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Figure 1. 

The eigenvalues E(X) of the Hamiltonian H = —A — 1/r + Xr for N = 3, n = 1, 
and £ = L = 0,1,2,3. The continuous curves show the upper bound EUL given 
by the envelope formula (2.9) with v = /i = P-j^(l), and the lower bound ELS 
by the sum approximation given by the same formula but with v = -Pj^(l) and 
[i = Pfi(— 1). The dashed curve EX represents accurate numerical data. 
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V(r)=-a/r+br. 



• h(r)= -a+p r 



V(r) 




Figure 2. 

The linear potential h(r) = ar + [3 used to estimate an upper bound for the 
eigenvalues of the Coulomb-plus-linear potential V(r) = —a/r + br . A and B are 
the absolute values of the inter-potential areas (or ip -weighted areas). We vary a 
and P so that A = B , and thereafter Theorems 5 and 6 imply E v < E h . 
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h(r)=a-p/r •* 


i i 


- 








— - V(r)=-a/r+br 
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*-A 
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i i 



12 3 4 
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Figure 3. 

The hydrogenic potential h(r) = —a/r + (3 used to estimate a lower bound for the 
eigenvalues of the Coulomb-plus-linear potential V(r) = —a/r + br . A and B are 
the absolute values of the inter-potential areas (or ifj -weighted areas). We vary a 
and P so that A = B , and thereafter Theorems 5 and 6 imply E h < E v . 
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Figure 4. 

We compare the bounds for E(X) , where E(X) is the ground-state eigenvalue 
(n = 1, I = 0) of the Hamiltonian H = —A — 1/r + Xr . The upper bounds 
(full-line) are by harmonic-oscillator tangents EUHO, linear tangents EUL, linear 
chords EUC, and linear chords with the wave function EUCW. The lower bounds 
(lower full-lines) are by hydrogenic tangents ELHY, Hydrogenic chords ELC, and 
Hydrogenic chords with the wave function ELCW. The dashed curve ELS represent 
the lower bound given by the sum approximation. Accurate numerical data (dotted- 
curve) EX is shown for comparison. 
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1 2 3 4 5 

X 

Figure 5. 

Bounds on the eigenvalues E^ (X) corresponding to the Coulomb-plus-linear po- 
tential V(r) = —1/r + Xr in N dimensions. Upper bounds EUC by the generalized 
comparison theorem, lower bounds ELS by the sum approximation, and accurate 
numerical data (dashed-line), for n = 1, £ = 0, and N = 3, 4, ... 7. By Theorem 2 
we know that the same curves apply also to i > since E^ e = E^ +2£ . 



